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A theoretical study of the intense-field single ionization of molecular hydrogen or deuterium ori- 
ented either parallel or perpendicular to a linear polarized laser pulse (400 nm) is performed for 
different internuclear separations and pulse lengths in an intensity range of (2 — 13) x 10^'^ Wcm~^. 
The investigation is based on a non-perturbative treatment that solves the full time-dependent 
Schrodinger equation of both correlated electrons within the fixed-nuclei and the dipole approxi- 
mation. The results for various internuclear separations are used to obtain the ionization yields of 
molecular hydrogen and deuterium in their ground vibrational states. An atomic model is used to 
identify the influence of the intrinsic diatomic two-center character of the problem. 



PACS numbers: 32.80.Rm, 33.80.Rv 



I. INTRODUCTION 

The theoretical and experimental investigation of the 
interaction of molecules with intense laser fields remains 
one of the most challenging problems of atomic, molecu- 
lar, and optical physics. Despite the numerous exper- 
imental work on molecules in intense ultrashort laser 
pulses during the last decades (see, e.g., P, 0| for re- 
views) a full understanding of the influence of the molec- 
ular structure on the strong-field response is still lack- 
ing. However, such an understanding is a prerequi- 
site for the recently proposed techniques that aim for 
the time-resolved imaging of changes of the electronic 
structure during chemical reactions. The development 
of such imaging techniques is driven by successful pi- 
oneering experiments in which the strong-field induced 
high-harmonic radiation Q or the in the ionization pro- 
cess ejected electrons Q were used to image the highest- 
occupied molecular orbitals of molecular nitrogen or oxy- 
gen. 

In order to achieve a three-dimensional image of the 
electron density or even an electronic orbital it is of 
course important to perform a spatially resolved mea- 
surement. Often, field-free alignment by means of rota- 
tional wave-packets [5, 6] is adopted for providing angle- 
resolved molecular strong-field data 0, 0, H Hi j but in 
[lol . [ll| alternative techniques without alignment are 
used. These experimental efforts are accompanied by 
a number of theoretical investigations of the orienta- 
tional dependent ionization probability of molecules in in- 
tense laser pulses. Since most molecular systems require 
some approximative treatment like the strong-field ap- 
proximation based additionally on the single- active elec- 
tron approximation (SAE) [l^j [3 [IB effec- 
tive independent-particle model like the time-dependent 
density- functional theory (TD-DFT) [l^, the simplest 
neutral stable molecule, H2 whose orientational depen- 
dence fiol [Til or angular distribution of ejected elec- 
trons [l8| was recently investigated experimentally is, in 
principle, an attractive alternative. Solutions of the time- 
dependent Schrodinger equation (TDSE) describing both 



electrons of H2 exposed to an intense la ser p ulse in full di- 
mensionality became recently available O |23, IIH . How- 
ever, these calculations were restricted to a parallel ori- 
entation of the molecule with respect to the field axis 
of a linear-polarized laser. This simplifies the treatment 
drastically, since the problem reduces to five spatial di- 
mensions, as the cylindrical symmetry is preserved. 

Theoretical investigations of the orientational depen- 
dence of the strong-field behavior of H2 are thus rather 
limited so far. This includes studies within the lowest- 
order perturbation theory (LOPT) [111, TD-DFT [H, 
or a Hartree-Fock based SAE approach [23] ■ The valid- 
ity of the SAE (and simplified models like the molecular 
Ammosov-Delone-Krainov tunneling model (MO-ADK) 
or the molecular strong-field approximation (MO-SFA)) 
that reduces the problem to three spatial dimensions was 
investigated in [2^ and found to be problematic espe- 
cially for few-photon processes. Only very recently the 
first investigation of the orientational dependence of the 
behavior of H2 in ultrashort intense laser pulses based 
on the TDSE was presented [265, and thus a full treat- 
ment of all six spatial dimensions of the two electrons 
was achieved. In [26] the single-ionization yield for a 
parallel and a perpendicular orientation was compared 
as a function of the wavelength, spanning an interval be- 
tween about 50 and 400 nm. Considering two different 
laser intensities and nuclear distances, only a brief idea 
of the influence of these parameters was provided. 

In the context of imaging it is, however, important to 
investigate in which way laser intensity or the quantum- 
mechanically unavoidable zero-point vibrational motion 
may blur the obtained image. On the other hand, 
the strong-field response itself may be used to visual- 
ize nuclear dynamics with sub-femtosecond resolution 
[13, m, [M [13, [m, m [33. In fact, for molecules like 
H2 a strong influence of nuclear motion on the strong- 
field ionization behavior was predicted on the basis of 
a simple model in [33] and later confirmed by ab initio 
calculations of quasistatic rates [35, 36, 37, 38] and full 
TDSE calculations [39| . An experimental confirmation 
for 800 nm radiation was achieved by the observation of 
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strong deviations from the Franck-Condon distribution 
of the formed vibrational states |40!] and the occur- 
rence of vibrational wavepackets in neutral H2 due to a 
phenomenon called Lochfraji [sTl that had been theoret- 
ically predicted in 30,] . However, the responsible strong 
dependence of the ionization yield on the intcrnuclear dis- 
tance was predicted for the so-called quasi-static regime, 
i. e. for low frequencies (long wavelengths) and high in- 
tensities. On the other hand, the perturbative results in 
[2^ indicated that the dependence on the internuclear 
distance is expected to be rather small in this so-called 
multiphoton regime (high frequency and low intensity). 

This work investigates the intensity and internuclear- 
distance dependence of the single-ionization yield of par- 
allel or perpendicular aligned H2 in ultrashort linear- 
polarized laser pulses with a wavelength of about 400 nm, 
as they are, e.g., experimentally available from a 
frequency doubling (second-harmonic generation) of a 
titanium-sapphire laser source. Besides the experimental 
relevance, the chosen wavelength is also of theoretical in- 
terest, since one expects six-photon ionization processes 
to dominate which lie somehow in the middle between 
few-photon and many-photon regimes. Thus one expects 
neither a simple perturbative nor the quasistatic approx- 
imation to be applicable. This ambivalent character is 
shown to be clearly visible, since, e.g., the ionization 
yield shows a pronounced dependence on the internuclear 
separation (as is expected for the quasistatic regime) , but 
also clear structures due to resonance-enhanced multi- 
photon ionization (REMPI). However, these structures 
are to a large extent washed out, if one goes beyond 
the fixed-nuclei approximation. In the present work vi- 
brational motion is considered in an approximate way 
where the internuclear-separation dependent ionization 
yields are weighted by the vibrational wavepacket of the 
initial state. The resulting ionization yields of H2 and 
D2 are compared with each other to resolve the isotope 
effect. Since for very short laser pulses the pulse dura- 
tion influences the strong-field behavior, the effects of the 
laser-pulse duration are investigated considering pulse 
lengths of 5, 10, and 20 fs. In the following atomic units 
(e — rrie = ft = 1) are used unless specified otherwise. 



II. METHOD 

Our method of solving the TDSE describing molecular 
hydrogen exposed to a laser field for parallel orientation 
and its extension to the perpendicular one is discussed in 
detail in [lO] and [111, respectively. Briefly, the TDSE is 
solved by expanding the time-dependent wave function 
in terms of field-free states. The latter are obtained from 
a configuration-interaction (CI) calculation [4l| in which 
the Slater determinants are formed with the aid of H^ 
wave functions expressed in terms of B splines in prolate 
spheroidal coordinates (1 < ^ < 00,— 1 < ?y < 1,0 < 
(j) < 27r). The use of a B-spline basis confined within a 
finite spatial volume defined by parameter ^max results 



in a suitable discretization of the electronic continuum. 



A. Configuration-interaction calculation 

For the considered laser parameters the electronic field- 
free states for every molecular symmetry were obtained 
as follows. A box size of about 350 oq is chosen indepen- 
dently of R. This is achieved by a proper adaption of ^max 
as a function of R. Along the ^ coordinate 350 B splines 
of order A; = 10 with an almost linear knot sequence were 
used where the first 40 intervals are increased by a ge- 
ometric progression using a progression factor g = 1.05 
and then the size of the interval remains constant. Along 
the 77 coordinate 30 B splines of order 8 were used in the 
complete interval —1 < ?7 < 1, but using the symme- 
try of a homonuclear system as is described in ,41] . Out 
of the resulting 5235 orbitals for every symmetry only 
3490 orbitals were further used to construct CI configu- 
rations, whereas those orbitals with highly oscillating an- 
gular part (with more than 19 nodes for the 77-dependent 
component) were omitted. In most of the subsequent CI 
calculations approximately 6000 configurations were used 
for every symmetry. These states result from a very long 
configuration series (3490 configurations) in which one 
electron occupies the Hj ground-state 1 ag orbital while 
the other one is occupying one of the remaining, e.g., 
nTTu or n Sg orbitals. The other CI configurations rep- 
resent doubly excited situations and are responsible for 
describing correlation (and real doubly excited states). 
Finally, out of the obtained CI states only those with an 
energy below the energy cut-off (chosen at 10 a.u. above 
the ionization threshold for the calculations shown in this 
work) were included in the time propagation (about 5400 
states per symmetry). For the perpendicular orientation 
only molecular symmetries with the absolute value of the 
component of the total angular momentum along the in- 
ternuclear axis < A < 7 were included in the time 
propagation. This results in a system of about 86,000 
real- valued first-order differential equations. Noteworthy, 
the adopted range of A values does not only guarantee 
the convergence of ionization yields, but also provides a 
reasonable description of photoelectron energy spectra. 

The basis set specified above was chosen to provide 
a good compromise for describing a large number of 
states and can, of course, not compete with a high- 
precision calculation optimized for a single electronic 
state. Fig. [1] demonstrates a comparison of the obtained 
electronic energies for different low-lying molecular states 
with high-precision calculations of L. Wolniewicz and 
co-workers performed using an explicitly correlated ba- 
sis 0, m, Si m, m. (The present CI method is able 
to reproduce such practically exact electronic energies at 
least within 4-6 significant digits, if the basis set is chosen 
judiciously [4l|). The agreement is excellent for all states 
except for l^S^ and 1^S+ (at larger internuclear dis- 
tances) where the electronic motion is highly correlated 
and cannot efficiently be described by a CI calculation 
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FIG. 1: Electronic potential curves of some low-l ying states of H2: the CI results of this work (blue circles) are compared to the 
supposedly exact values (solid lines, [i^. |45|. |46| ) . The dashed line shows the ionization threshold of H2. The different 

adiabatic electronic states of a given molecular symmetry are numbered in the order of the energy. 



employing orbitals with no electron-electron interaction 
included. Nevertheless, even for these two states the ob- 
tained electronic energies are much better than those ob- 
tained with the Hartree-Fock approximation. For exam- 
ple, for the ground state of H2 with the exact electronic 
energy at i? = 1.4 ao being equal to —1.1745 a.u., the 
Hartree-Fock limit is —1.1336 a.u., whereas the present 
CI calculation yields —1.1604 a.u. 



tential curve due to the external field. The latter assump- 
tion implies that the index v of vibrational state should 
be sufficiently small. In this work, only = is con- 
sidered, but for the two isotopes H2 and D2 with their 
different vibrational wavefunctions. Due to the larger 
mass the D2 vibrational ground state is more compact 
than the one of H2 . 



B. Integration over internuclear separations 

Once the TDSE is solved for a given linear-polarized 
laser pulse, fixed internuclear separation i?, and angle d 
between the internuclear and the polarization axis, the 
ionization yield Y\a^^(I{^ 0) is obtained from a summation 
over the populations of all discretized continuum states. 
These yields can be further used to calculate the ion- 
ization yield i^o^H^) ^ given initial vibrational state 
v described by the vibrational wavefunction cfy^i^R). In- 
deed, if the duration of the pulse is sufficiently short and 
depletion of the state during the pulse is insignificant, one 
can neglect the motion of the wavepacket created during 
the pulse and calculate the total ionization yield as [s^l 



di?yio„(i?, e)\(j),{R)[' 



(1) 



Application of Eq. ^ is further based on the assumption 
that the molecule has no time to rotate during the pulse 
and neglect distortion of the electronic ground-state po- 



C. Atomic model 

For the analysis of the orientational dependence of the 
ionization due to the anisotropy of a molecule, it is con- 
venient to compare the molecular results with those ob- 
tained for an artificial atom with an isotropic, single- 
centered charge distribution. Strong-field ionization is, 
however, known to be not only sensitive to the symme- 
try, but also to the electronic binding energy and the 
exact form of the long-ranged Coulomb potential. There- 
fore, the artificial atom must agree to the corresponding 
molecule with respect to the two latter factors. For this 
reason the simple single-electron one-parameter model 
potential 



V{r) 
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1 
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2r 



11/2 



(2) 



was introduced [26|. Its performance for describing var- 
ious physical problems was checked in [i^ and it was 
recently also applied to the calculation of antiproton-H2 
scattering cross sections and stopping powers [i^, l49j . 
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FIG. 2: Dependence of the atomic-motiel parameter a [see 
Eq. ^] on the intemuclear distance R. This choice of a(7?) 
leads to a vertical ionization potential of H2 in agreement with 
the one of the present CI calculation. 



The ionization potential Ip of such an artificial atom is 
directly related to the parameter a of the model potential 
(|2l). For a given ionization potential Ip the correspond- 
ing parameter a can be found as a = a{Ip), where a{Ip) 
can be obtained numerically. Since the molecular ver- 
tical ionization potential (energy difference between the 
electronic ground-state potential curves of the ion and 
the neutral) depends on the intemuclear distance R, the 
value of a should also depend on R in order to compare 
atomic- model and molecular results. Figure [2] shows the 
R dependence of a. It is determined by requiring the 
resulting i?-dependent ionization potential to agree with 
the one obtained by the present CI calculation for the II2 
ground state. 

In order to compare to the full molecular calculations 
the atomic-model results are multiplied by a factor 2 that 
accounts for the two equivalent electrons in molecular 
hydrogen. This procedure is known to be reasonable for 
ionization yields less than 10-20% [2^ which is the case 
for the present calculations. 



III. RESULTS 

All calculations presented in this work were performed 
with A^-cycle cos^-shaped linear-polarized laser pulses 
with N = 10,20, and 40. For a wavelength of 400 nm 
the FWHM of intensity of such pulses corresponds to 
about 5, 10, and 20 fs. 



A. Field-induced resonances 

For the following discussion of the orientation depen- 
dence and isotope effects, it is helpful to obtain a more 
detailed understanding of the influence of such param- 
eters as peak intensity and wavelength of the pulse or 



FIG. 3: Expected positions of REMPI peaks and A'^-photon 
ionization thresholds. Dependence of the resonant intensities 
/ (a) on the intemuclear distance R for pulses with a wave- 
length of 400 nm and (b) on wavelength for the fixed intemu- 
clear distance R = 1.4 ao. Positions of REMPI peaks due to 
(red solid), n„ (blue dashes), or (black solid) inter- 
mediate states are shown. Also shown are the Af-photon ion- 
ization thresholds (dotted lines, with A'' explicitly indicated 
in the graph). 



of the intemuclear separation on the positions of REMPI 
peaks and TV-photon ionization thresholds. In order to be 
able to correctly predict REMPI through some resonant 
electronic state in intense laser pulses, it is necessary to 
know the field-induced shift of the resonant state, what 
is a challenging task by itself. However, if the field is suf- 
ficiently intense, one can assume that the field-induced 
shift (dynamically induced Stark shift) of excited states 
is almost equal to the ponderomotive energy. With this 
assumption and using the field-free transition energies of 
the present CI calculation the positions of the REMPI 
peaks are expected to depend on the laser parameters 
and intemuclear separation as shown in Fig. [3l 

The dependence of the expected positions of the 
REMPI peaks and iV-photon ionization thresholds on 
the intemuclear distance R for a 400 nm laser field is 
given in Fig. [3]a. In this R range the ionization process 
can be referred to as 5-photon (7-photon) ionization in 
the bottom-right (top-left) part of the figure, or as 6- 
photon ionization otherwise. Different kinds of REMPI 
peaks are expected: (5-1-1) REMPI peaks through nS+ or 
tiHu electronic states with n > 1, (5+2) REMPI peaks 
through the 11]+ or in„ states, a (4+1) REMPI peak 
through the 3S+ state, and a REMPI peak through the 
2S+ state. Note, that in the last case the expected posi- 
tion of the resonance crosses the expected position of the 
5-photon ionization threshold. Therefore, the resonance 
can be referred to as (4+1) REMPI for peak intensities 
smaller than 7.5 x 10" W/cm^, and as (4+2) REMPI for 
higher intensities. Evidently, the correct character of the 
resonance is sensitive to the exact intensity dependence 
of the field-induced shift of the 2E+ state, and thus a 
non-trivial behavior is expected. A similar conclusion is 
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(a) i?, = 1.4ao (b) R = 2.2ao 

FIG. 4: Wavelength-dependent ionization yields for a parallel 
orientation of an H2 molecule at the fixed internuclear dis- 
tances (a) R = 1.40 ao and (b) K = 2.20 ao for various peak 
intensities (specified in units of 10^'^ W/cm'^) of 40-cycle cos^- 
shaped pulses. The arrows indicate the resonant wavelengths 
(in the low-intensity limit) of REMPI peaks due to nE^^ (a) 
and nE^ (b) intermediate states. The expected (dashes) and 
found (dots) position of a REMPI peak is also given in (a). 

valid for the REMPI through the 1E+ or in„ electronic 
states, since their exact REMPI positions could in fact 
cross the 6-photon ionization threshold. 

Similarly, Fig. Ob shows the dependence of expected 
positions of REMPI peaks and TV-photon ionization 
thresholds on the laser wavelength for the fixed internu- 
clear distance i? = 1.4 ao. With larger peak intensity the 
increasing ponderomotive energy leads to an increase of 
the transition energy between the initial and the resonant 
state. This increase can be compensated by the increase 
of the photon energy, and thus the new REMPI position 
will occur at a smaller wavelength. The calculations at 
a fixed internuclear separation are more suitable for the 
investigation of the validity of the assumed field-induced 
shift of electronic states, since in this case identical sets 
of field-free electronic wavefunctions are used in the time 
propagation. For this purpose, a series of 210 calculations 
for a parallel-oriented H2 molecule with fixed internuclear 
distance R = 1.4 ao exposed to laser pulses with a total 
duration of 40 cycles was performed for 21 different val- 
ues of the wavelength and 10 different values of the peak 
intensity. The results are shown in Fig. |3]a, where every 
point represents the outcome of one full TDSE calcula- 
tion, and curves join the results obtained for the same 
peak intensity. 

Figure [4]a shows a pronounced peak whose position 
moves from 407 nm for a peak intensity of 5 x 10^^ W/cm^ 
to 387 nm for a peak intensity of 7 x 10^^ W/cm^. From 
Fig. Ob it follows that this peak can be assigned to 
REMPI through either one or both of the closely ly- 
ing 3E+ and 4E+ electronic states, since the spectral 
width of the Fourier-limited pulse is too broad to resolve 
these two resonances. Clearly, the intensity-dependent 



shift of the peak position is overestimated by the already 
mentioned simple prediction based purely on the pon- 
deromotive energy {SE{I) = Up). Instead, the found in- 
tensity dependence of the field-induced energy shift can 
be well fitted by SE{I) = 0.9Up - 0.002. At intensi- 
ties 5 X 10^^ W/cm^ and smaller it appears as the energy 
shift of these low-lying excited states (responsible for the 
REMPI) is aheady absent. The position of the REMPI 
peak agrees then much better with the low-intensity limit 
than with the prediction based on Up, since the latter 
would predict a shift of about 2 nm. Such a shift by 
about 2 nm is, however, found for the (poorly resolved) 
REMPI peak due to the higher lying intermediate 
state. This demonstrates that in the investigated regime 
of laser parameters different excited states behave differ- 
ently, and a common prediction for all excited states is 
impossible. 

Figure |3]b shows again the results of a series of (this 
time 120) calculations for a parallel-oriented II2 molecule, 
but for the larger internuclear separation R = 2.2 ao. In 
this case the spacing of the 3E+ and 4E+ states that 
could lead to (4-1-1) REMPI is rather large, and thus 
the pronounced peak in Fig. [4]b can be entirely assigned 
to REMPI through the 3S+ state. Although the peak 
position clearly shifts to smaller wavelengths with in- 
creasing laser peak intensity, the shift becomes visible 
only for rather large intensities. As a consequence, the 
peak position crosses 400 nm at an intensity higher than 
10^"^ W/cm^, whereas according to Fig. [3]a the crossing 
should have occurred at an intensity that is smaller by a 
factor 2. Thus, although the positions presented in Fig. [3] 
give a satisfactory explanation of the main features, they 
should only be considered as a rough estimate. 



B. J?-dependent ionization 

If strong-field ionization of II2 or D2 initially in their 
vibrational ground states is considered, it is important 
to investigate the dependence of the ionization on the 
internuclear separation R within an R range in which 
the vibrational wave function is nonvanishing (Franck- 
Condon window). Therefore, the TDSE describing II2 
within the fixed-nuclei approximation was solved for 25 
different values of R (in between 1.0 ao and 2.2 ao with a 
step size of 0.05 ao). This was repeated for 15 different 
values of peak intensities (in a range from 2 x lO^'^W/cm^ 
to 1.3 X lO^-^W/cm^). 

Figure [5]a (left panel) shows the obtained results for a 
parallel orientation of the molecular axis with respect to 
the polarization vector and 40-cycle cos^-shaped (FWHM 
of about 20 fs) laser pulses. Two main features may be 
observed. First, the ionization yield increases with R. 
Second, pronounced structures are visible. The impor- 
tance of both effects decreases with intensity. The in- 
crease with R was first predicted in the quasistatic regime 
[3^ 1. Its occurrence at 400 nm shows that even for this 
wavelength clear strong-field phenomena are observable 
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(a)Parallel orientation of the internuclear axis to a linear polarized laser field 




Intemuclear distance R (a.u.) Internucleai' distance R (a.u.) Internuclear distance R (a.u.) 



(b) Perpendicular orientation of the internuclear axis to a linear polarized laser field 




FIG. 5: Ionization yields for 40-cycle cos^-shaped (20 fs) laser pulses with a wavelength of 400 nm and different peak intensities 
in between 2.0 and 13 xlO^^ W/cm^ (specified in the graphs in units of 10^'^ W/cm^). The upper (middle) panel shows the 
results for a parallel (perpendicular) orientation of the molecule with respect to the field axis, while the lower panel shows 
the results obtained with the atomic model potential in Eq. ([2]). The left panel shows the fixed-nuclei ionization yields as a 
function of the internuclear distance R, whereas the middle and right panels display the ionization yields multiplied with the 
probability density of the ground vibrational states of H2 and D2, respectively. (Every point corresponds to a full solution of 
the TDSE, while the points are connected by spline interpolating curves.) 
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(a)Parallel orientation of the intemuclear axis to a linear polarized laser field 
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(b) Perpendicular orientation of the intemuclear axis to a linear polarized laser field 
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(c)Atomic model calculations (multiplied by a factor 2) 
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FIG. 6: As Figure O but for a 20-cycle cos^-shaped (10 fs) laser pulse. 
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(a)Parallel orientation of the intemuclear axis to a linear polarized laser field 
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(b) Perpendicular orientation of the intemuclear axis to a linear polarized laser field 
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(c)Atomic model calculations (multiplied by a factor 2) 



FIG. 7: As Figure [S] but for a 10-cycle cos^-shaped (5fs) laser pulse. 



beyond a pure multiphoton picture. On the other hand, 
the observed structures are due to classical multiphoton 
phenomena (channel closings and REMPI). 

Some pronounced REMPI peaks are visible that should 
be compared with their predicted positions in Fig. [31 The 



peaks at i? = 1.2 oq and i? = 1.4 oq for a laser peak in- 
tensity of 2 X lO^^W/cm^ can thus be assigned to (5-1-1) 
REMPI through the 2 E+ and 3 - 4 E+ states, respec- 
tively. The position of the latter peak changes with 
the peak intensity almost in the expected way as has 
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also been demonstrated in Fig. |4]a. On the other hand, 
the position and amphtude of the REMPI peak arising 
from the 2 resonant state cannot easily be under- 
stood. For a peak intensity of lO^^'W/cm^ this REMPI 
peak is located at R = 1.6 oq instead of the expected 
value of 1.5 ao- Interestingly, the amplitude of the peak 
that is very large at small intensities is becoming very 
small for higher intensities, as one may expect when go- 
ing from the multiphoton in the direction of the qua- 
sistatic regime. For intensities between 4 x 10^'^W/cm^ 
and 6 x lO^^W/cm^ one can observe something similar 
to a splitting of the REMPI peak into two peaks. This 
behavior can evidently not be explained using Fig. [S] 
It is also difficult to explain the pronounced peak lo- 
cated at R — 1.4 ao for the highest laser peak intensity 
(1.3 X 10^'^W/cm^). According to Fig.Ono peak should 
occur for these values of R and /, since it lies in between 
the expected positions of the 1 E+ and the 2 S+ REMPI 
peaks. This may indicate some field-induced coupling of 
these states and thus a clear strong-field phenomenon. It 
is also interesting to note that the channel closings indi- 
cating the transitions from 5- to 6-photon ionization and 
from 6- to 7-photon ionization are visible, but not very 
pronounced. Furthermore, due to the REMPI peaks the 
channel thresholds are sometimes difficult to identify in 
the shown ionization yields. 

In order to consider the influence of vibrational mo- 
tion onto the strong-field ionization yields, the results ob- 
tained for a fixed nuclear orientation are weighted with 
the probability density of the ground vibrational state 
[see integrand of Eq. The corresponding result for 

H2 is shown in the middle panel of Fig.O Only in the case 
of the highest laser peak intensity considered in this work, 
the maximum of the weighted ionization yield agrees 
with the maximum of the vibrational wave function. At 
slightly lower intensities (until about 8 x 10^'^W/cm^) 
the weighted ionization yield is largest for smaller values 
of R, while for even lower intensities the REMPI peaks 
due to the 3S+ and 4E+ states determine the maxi- 
mum of the weighted ionization yield. At the lowest laser 
peak intensity considered (2 x 10^'^W/cm^) the highest 
weighted ionization yield is found at around R = 1.7 oq. 
This maximum should be due to an opening of the 5- 
photon regime and may be further increased by (4-1-1) 
REMPI processes. At this intensity one notices also a 
very slow decrease of the weighted ionization yield for R 
values above 2ao, despite the fact that the vibrational 
wave function has a very small amplitude. The reason is 
the aheady discussed 3 E+ REMPI peak (Fig.[4]b). 

In the case of D2 (left panel of Fig. [5]a) the narrower 
vibrational distribution is, however, sufficient to domi- 
nate over the resonant effect. In this case the weighted 
ionization yield decreases rather pronouncedly for inter- 
nuclear separations larger than 2.0 ao- Close to the min- 
imum of the electronic potential curve at about 1.4 
the weighted ionization yields are, however, very simi- 
lar for H2 and D2 . Also for D2 the ionization yield peaks 
only for the highest intensity considered here at the max- 



imum of the vibrational wavefunction. Consequently, the 
Franck-Condon approximation would not describe the vi- 
brational distribution of the formed ions properly. A 
proper calculation of these distributions has to include 
the effects of channel closings and REMPI, but also of 
the general increase of the ion yield as a function of in- 
ternuclear separation. Neither a pure multiphoton nor 
quasistatic prediction is thus sufficient. 



C. Orientational dependence 

For a perpendicular orientation of the molecule with 
respect to the field (Fig. Ob) the i?-dependent ionization 
yield (left panel) looks on the first glance surprisingly 
similar to the one for the parallel orientation. For the 
lowest shown intensities the spectra comprise very pro- 
nounced peaks at about R = 1.25 ao and R — 1.4 ao 
that can be assigned to (5-1-1) REMPI through the 2n„ 
and the 3 — 4n„ states, respectively. It is a peculiar- 
ity of II2 that already the lowest lying excited states of 
and ^IIu symmetry and thus REMPI peaks through 
those states lie energetically very close together. In con- 
trast to the results for parallel orientation a third peak 
at i? = 1.55 oq is, however, also well resolved. From 
Fig. [3] it appears very likely that this peak stems from 
a superposition of (5-1-1) REMPI processes through the 
higher excited n„ states. Especially at lower intensities 
one notices furthermore that the ionization yield does not 
increase that evidently for large R values than it does for 
a parallel orientation. The threshold between 5- and 6- 
photon ionization is rather well resolved and appears for 
the different laser peak intensities more or less at the 
expected R values (Fig. [3]a). As a consequence of the 
smaller slope at large R the weighted ionization yield 
decays for larger R values much faster for the perpendic- 
ular than for the parallel orientation for both H2 (middle 
panel of Fig. [Sib) and D2 (right panel). 

Figure [5]c shows finally the ionization yield obtained 
with the simple isotropic one-electron model potential 
given in Eq. The agreement of the i?-dependent 

ionization yields obtained with this model and the full 
molecular two-electron calculation are surprisingly good, 
especially with the results obtained for the parallel ori- 
entation. As in the latter case, the atomic model gives 
a shifted threshold between the 5-and 6-photon regimes 
compared to the prediction according to Fig. [3l The 
atomic model yields also a rather pronounced increase 
in ionization for large R values, especially for low laser 
peak intensities as was also found for the parallel ori- 
entation. A closer look reveals, however, that for small 
R separations and especially for the first REMPI peak 
the atomic model agrees slightly better with the molec- 
ular results obtained for the perpendicular orientation. 
The 2nd REMPI peak defines somehow the transition 
line. For smaller R values the atomic model agrees bet- 
ter with the perpendicular results, while starting with the 
2nd REMPI peak the ionization yields obtained for the 
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FIG. 8: Comparison of ionization yields (in a 20 fs laser pulse) for different internuclear distances R (weighted with the the 
probability density of the ground vibrational state of H2 at this R) for a parallel (red curves) and a perpendicular (blue curves) 
oriented H2 molecule. The R values are specified explicitly in the figure. The vertical dashed lines indicate the position of a 
peak intensity 6 x lO'^^W/cm'^ on the intensity axis for every corresponding pair of curves. 



atomic model and the molecular one for a parallel orien- 
tation agree better with each other. The main difference 
to the molecular calculations is the position of the 3rd 
REMPI peak that for the lowest shown intensity lies so 
close to the 2nd one, that it appears in the i?-dependent 
ionization yield as a shoulder. For a laser peak inten- 
sity of 2.5 X lO^^-y^/cm^ the 3rd REMPI peak is shifted 
more than the 2nd one and is thus visible as a well sepa- 
rated peak. However, for higher intensities it is less well 
resolved due to its low probability. Despite the overall 
good agreement of the results for the atomic model with 
the full molecular calculations (on a logarithmic scale!), 
the weighted ionization yields still reveal differences. For 
example, the maximum of the weighted ionization yields 
for H2 and D2 and the largest laser peak intensities is 
shifted to slightly smaller R values than is found for the 
full molecular calculations. 

A further important laser parameter is the pulse dura- 
tion. Its influence is demonstrated in Figs. [5] and [7] that 
show the corresponding results for 20- and 10-cycle pulses 
(FWHM of 10 and 5fs), respectively. The increased laser 
bandwidth leads to spectra that show much less details 
compared to the relatively long 40-cycle pulse. The ion- 
ization yields for the 10-cycle pulse show almost no evi- 
dence of REMPI peaks. The curves are fairly smooth and 
the remaining structures can be explained by the closing 
and opening of iV-photon ionization channels. 

Because of the different positions of the REMPI peaks 
for parallel or perpendicular orientations the ratio of 
parallel to perpendicular ionization yields may substan- 
tially change for a small variation of R. This effect is 
demonstrated in Fig. [8] in which ionization yields (mul- 
tiplied with the probability density of the ground vi- 
brational state) for parallel and perpendicular oriented 
H2 molecule are compared for different internuclear dis- 



tances. A log-log scale is used and the pairs of curves 
(parallel and perpendicular orientation for a given value 
of R) are shifted along the intensity axis for better read- 
ability. To guide the eye, the vertical dashed lines indi- 
cate the position of the peak intensity 6 x 10^'^W/cm^ 
on the intensity axis for every pair of curves. The multi- 
plication with the probability density was performed in 
order to emphasize the relative contributions of different 
internuclear distances for the total ionization ratio be- 
tween parallel and perpendicular orientations obtained 
after integration over R. 

As can be seen from Fig. [51 the ionization yields for 
parallel and perpendicular orientation are almost equal 
in the range R = 1.3 — 1.7 oq for a peak intensity of 
6 X lO^^W/cm^. Whereas the ionization yield for the 
parallel orientation is larger than for the perpendicu- 
lar one at i? = 1.3 and 1.5 ag, the opposite is found at 
R = 1.4, 1.45, and 1.6 oq. At smaller values of R the par- 
allel orientation is slightly easier ionized than the perpen- 
dicular one, while for R values larger than 1.8 oq parallel 
oriented molecules are much easier ionized. This is a 
consequence of the slower decay of the ionization yield 
for a parallel orientation and for large R values that was 
already discussed in the context of Figs. [5]a and b. 

The key conclusion that can be drawn from Fig. [H] 
is the need for systematic studies of the intensity and 
internuclear-separation dependencies of the ratio be- 
tween the ionization yields for parallel or perpendicular 
orientation as they are performed in this work, since a 
calculation for a single laser peak intensity and internu- 
clear separation R can yield any possible result, i. e. the 
ratio between the ionization yields for parallel and per- 
pendicular orientations may be found to be equal to 1, 
much smaller than 1, or much larger than 1. Depending 
on the choice of intensity and R very different conclusions 
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on the orientation dependence of the ionization yield of 
H2 in strong laser fields would follow. 

D. Integrated ionization yields 

The i?-independent ionization yields for H2 and D2 
molecules (in their vibrational ground states) are finally 
obtained by an integration of the weighted i?-dependent 
ionization yields over R (Eq. ^ in Sec. IIIB|) . The re- 
sults for parallel and perpendicular orientation are shown 
for the diS'erent pulse lengths in Fig. [SI The curves ob- 
tained after R integration are much less structured than 
the fixed-nuclei ionization yields, as is evident from a 
comparison with Fig. [8|). Clearly, the structures due to 
REMPI processes are smoothed out by the integration 
over R. As a consequence, the curves look almost like 
straight lines on the used log-log scale. 

In fact, it turns out that the linear dependence on the 
logarithmic scale is well described, if the yield is fitted 
with the function 

Y,{I,T) = nT{I/Io)''' (3) 

where / is the peak intensity, Iq = 3.5094452 x 10^^ 
W/cm^ is the atomic unit of intensity, and T is the 
FWHM duration of the pulse in atomic units. For the 
fit parameters the values = 1.55 ■ 10^ and kg = 4.17 
are found. The obtained value of kg indicates a non- 
perturbative behavior, since according to Fig. [3]a one 
would expect mostly 6-photon ionization to occur and 
thus ks should be close to 6. The included dependence 
on the pulse duration T allows to compare the results 
obtained for different pulse lengths. A linear dependence 
on T should be found, if a rate concept is applicable. 

Dividing the ionization yields by the fit function Q 
allows a direct comparison of the parallel, perpendicu- 
lar, and atomic model potential results on a linear scale. 
They are shown in the inserts of Fig. [9l If the vibra- 
tional ground state of II2 is considered (middle panel of 
Fig. [9|) , the scaled atomic-model results are closest to 1 
and thus most accurately described by the fit function. 
The scaled yield for parallel orientation decreases from 
a value of 1.5 to about 1.0 for the 10-cycle pulse, but 
shows an increasing behavior for higher intensities in the 
case of the 20-cycle pulse. The smallest intensity depen- 
dence is found for the longest pulse considered in this 
work where the scaled yield varies only between about 
1.25 and 1.15. Interestingly, the scaled yield for perpen- 
dicular orientation shows almost the opposite behavior. 
The most pronounced intensity dependence is found for 
the longest pulse. Furthermore, the scaled yield increases 
for low intensities as a function of intensity. As a conse- 
quence, the scaled yields for parallel and perpendicular 
orientation first approach each other before they separate 
again for even larger intensities. 

Using the same fit function for scaling the D2 yields 
one notices that the yields for parallel orientation are 
now almost flat (shortest pulse) or increase with intensity. 



Since the vibrational density distribution of D2 is more 
localized around Rq, one can conclude that such higher 
ratio for H2 can be due to contributions to the ionization 
from either small or large internuclear distances R. The 
comparison of the middle and the right panels of Fig. [7]a 
shows that the effect stems from the enhanced ionization 
at R > 1.6 flQ. Also the scaled yields for perpendicular 
orientation or the atomic model potential show a larger 
increase with intensity for D2 compared to H2 , although 
this effect is a little bit less pronounced. This indicates 
that the ionization yield of D2 possesses a slightly larger 
slope than the one of H2, a rather unexpected (though 
small) isotope effect. 

Finally, the ratio of the ionization yields for parallel 
to perpendicular orientation of the molecular axis as a 
function of the peak intensity is also shown in Fig. [51 
For the 10-cycle pulse and H2 the ratio is about 2 for the 
peak intensity 2 x 10^'^W/cm^ and decreases smoothly to 
about 1.1 at 9 X lO^^W/cm^, before it increases to 1.18 
at 1.3 X 10^*W/cm^. The occurrence of a minimum is 
due to the maximum found for scaled yield in the case 
of the perpendicular orientation, as was discussed in the 
context of the inserts in the middle and right panels of 
Fig. [21 Increasing the pulse duration docs not change the 
behavior in a qualitative fashion, but the ratio found at 
small intensities decreases with increasing pulse length. 
At the same time, the increase at the highest intensities 
is more pronounced, but this increase is smaller than the 
decrease seen for the low intensities. In the case of the 
40-cycle pulse the ratio starts at about 1.5, decreases 
to almost 1.0 and increases to 1.26. The turning point 
shifts also to slightly lower intensities for longer pulses. 
The intensity dependence of the ratios for D2 show a 
similar behavior as was found for II2. However, at small 
intensities the ratio is clearly smaller than for II2, i.e. 
the anisotropy of the ionization yield is less pronounced. 
For high intensities the ratios found for H2 and D2 agree 
on the other hand almost perfectly with each other. An 
isotope effect occurs thus only for low intensities. 

The importance of the inclusion of nuclear motion is 
evident from the ratio of parallel to perpendicular ion- 
ization yields obtained for a fixed internuclear separation 
{R = 1.4 ao) that is also shown in Fig. [9l In the case of 
a 10-cycle pulse the ratio is also smooth, but increases 
with intensity. The pronounced decrease found for the 
i?-integrated ratio at low intensities is thus completely 
absent. Interestingly, the agreement with the ratio found 
for II2 is very good for high intensities for which also the 
II2 and D2 ratios agreed well with each other. This indi- 
cates that for high intensities the ratio is less sensitive to 
R. The reason is the less pronounced R dependence of 
the ionization yields for high intensities that was found 
in general and discussed in the context of the weighted 
i?-dependent ionization yields in Figs. [5l to [3 One of the 
consequences of this reduced R dependence was, e.g., 
that the maximum of the ionization yield was more or 
less found for the R value at which the vibrational den- 
sity had its maximum. For longer pulses pronounced 
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FIG. 9: Final ionization yields (integrated over the internuclear separation R) as a function of the laser peak intensity for H2 
(middle panel) and D2 (right panel) and a parallel (red circles) or a perpendicular (blue squares) orientation and a) 10-cycle, b) 
20-cycle, or c) 40-cycle laser pulses. Also shown are the results obtained with the isotropic one-electron model potential (green 
triangles). The inserts show the ratio of the ionization yields to the fit function in Eq. (|3}- The resulting ratio of parallel to 
perpendicular ionization yields is shown in the left panel in which also the corresponding ratio obtained for a fixed internuclear 
separation {R — 1.4 ao) is plotted. 
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intensity-dependent maxima and minima become visible. 
As a consequence, the ratio found for a single R value 
differs clearly from the i?-integrated results. For exam- 
ple, in the case of a 40-cycle pulse and R = 1.4 ao the 
ratio decreases to about 0.5 at a laser peak intensity of 
5 X lO^^W/cm^ which means that perpendicular oriented 
H2 ionizes much better than parallel oriented one. This 
is in complete contrast to the iZ-integrated results for 
which the perpendicular orientation never ionizes faster 
than parallel oriented molecules in the considered range 
of laser peak intensities. 

IV. CONCLUSION 

We have performed an extensive numerical study of 
intense-field ionization of molecular hydrogen and deu- 
terium numerically integrating the full-dimensional two- 
electron Schrodinger equation in the non-relativistic, 
fixed-nuclei, and dipole approximation. The pre- 
sented results are obtained for three different durations 
(5, 10,20fs) of ultrashort frequency-doubled Ti:sapphire 
laser pulses (400 nm) for both parallel and perpendic- 
ular orientations of the molecular axis with respect to 
the laser field. Calculations are performed for 15 dif- 
ferent intensities (in a range from 2 x lO'^'^Wcm"^ to 
1.3 X 10^^ Wcm^^) and 25 different internuclear separa- 
tions (in a range from 1.0 ag to 2.2 ag) which results in 
375 data points for each orientation. The same series of 
calculations was performed employing an isotropic one- 
electron model potential in order to study the influence 
of molecular anisotropy or the one of the two electrons. 

By analyzing the dependence of the fixed-nuclei ion- 
ization yields on the peak intensity and the internuclear 
separation, we assign the observed peaks to REMPl or 
closings of A''-photon ionization channels and study field- 
induced shifts of resonant electronic states. 

A key feature of the present work is the calculation 
of total ionization yields of molecular hydrogen and deu- 
terium in their ground vibrational states by an integra- 
tion of the fixed-nuclei ionization yields multiplied with 
the corresponding vibrational probability density over 



the internuclear separation. The subsequent analysis of 
the ionization anisotropy reveals a smooth dependence 
on the laser peak intensity and pulse duration. The ob- 
tained ratios of the parallel to perpendicular total ion- 
ization yields vary in the investigated intensity and pulse 
duration ranges in between 1 and 2. Whereas for high 
intensities the ionization anisotropy of H2 and D2 is al- 
most identical, the ratio for H2 is always larger than the 
one for D2 in the case of small intensities, although the 
difference does not exceed 20%. 

The importance of the integration over the internu- 
clear separation for the interpretation of experiments is 
demonstrated by comparing the obtained parallel to per- 
pendicular ratios with those calculated using the fixed- 
nuclei results at the equilibrium distance. For example, 
there is a wide range of intensities where at the equilib- 
rium distance the ionization for perpendicular-oriented 
molecule is larger than for parallel-oriented ones. Be- 
sides the rather different qualitative dependence on the 
peak intensity, also the dependence on the pulse duration 
is different. 

The results of this work do not only provide an interest- 
ing insight into the strong-field behavior of molecules and 
the influence of channel closings, REMPl, field-induced 
shifts of energy positions, pulse length and intensity, vi- 
brational motion, isotope effects, and orientation, but 
will hopefully also stimulate experimental efforts to mea- 
sure the anisotropy of the ionization yield of H2 (or D2) 
in the interesting regime that is neither well described 
by the multiphoton nor the quasistatic approximations. 
Furthermore, the results should serve as benchmarks for 
other theoretical approaches and simplified models. 
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